NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


DTIC 

ELECTE 


JUN6  £66 


I 


EVALUATION  OF  THE  EFFECTS  OF  ROLL  ANGLE  ON 
THE  NATURAL  CONVECTION  HEAT  TRANSFER  FROM  AN 
ARRAY  OF  CYLINDRICAL  SPINES 

by 

Joseph  E.  McClanahan 
December  1984 

Thesis  Advisor:  A.  D.  Kraus 


Approved  for  public  release;  distribution  unlimited 


85  5  09  073 


security  classification  of  this  page  (TOim  dmi  Entered) 


REPORT  DOCUMENTATION  PAGE 


REPORT  number 


4.  title  (tnd  Submit)  _  '  „  ,  ,  , 

Evaluation  of  the  Effects  of  Roll  Angle 
on  the  Natural  Convection  Heat  Transfer 
from  an  Array  of  Cylindrical  Spines 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


«.  PERFORMING  ORG.  REPORT  NUMBER 


AUTHORS 

Joseph  E.  McClanahan 


i.  CONTRACT  OR  GRANT  NUMBERS A) 


>.  PERFORMING  ORGANIZATION  NAME  ANO  AODRESS 

Naval  Postgraduate  School 
Monterey,  California  93943 


to.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  *  WORK  UNIT  NUMBERS 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93943 


12.  REPORT  DATE 

December  1984 


13.  NUMBER  OF  PAGES  _ 


4.  MONITORING  AGENCY  NAME  6  ADDRESS  (it  dlttoront  from  Controlling  Oftlem)  IS.  SECURITY  CLASS,  (ot  thtm  roporl) 

Unclassified 


ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


•  6.  DISTRIBUTION  STATEMENT  (ot  thlt  1 Import ) 


Approved  for  public  release '  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (at  the  tbetrtct  mntmrmd  In  Block  20,  II  different  from  Rap 


.  O  CO 
!t»  I—  ^ 


IS-  KEY  WORDS  (Continue  on  rmvtrtm  tide  It  nacaa.ary  and  Identity  by  bloc*  number) 

/natural  convection',  J'inned  array ;  -cylindrical  spines,  heat  trans¬ 
fer,  roll  angle  enhancement Nusselt  number,  Rayleigh  number 


20  ABSTRACT  (Continue  on  r«v« r««  «fd«  If  n«c«*««ry  and  Idontity  by  block  miflibor) 

Experimental  data  has  been  obtained  that  predicts  the  effect  of 
roll  angle  upon  the  natural  convection  heat  transfer  from  an  ar¬ 
ray  of  cylindrical  spines  or  pin  fins.  A  correlation  is  present¬ 
ed  that  permits  evaluation  of  the  Nusselt  Modulus  as  a  function 
of  the  Rayleigh  Modulus  and  the  roll  angle  and  it  shows  that  max¬ 
imum  1  tat  transfer  is  obtained  at  roll  angles  between  20  and  30 
degrees . 


DO  /j 


FORM 
JAN  71 


EDITION  OF  1  NOV  41  IS  OBSOLETE 

S  N  0102-  LF-  014-  6601 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fBTiwi  Dot*  Bntortd) 


Approved  for  public  release;  distribution  is  unlimited 


Evaluation  of  the  Effects  of  Boll  Angle 
on  the  Vatural  Convection  Heat  Transfer  from  an  Array 

of  Cylindrical  Spines 


by 


Joseph  E.  McClanahan 

Lieutenant  Commander.  United  States  Navy 
B.S.,  Auburn  University,  1970 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
December  1984 


ABSTRACT 


~7  Experimental  data  has  been  obtained  that  predicts  the 
effect  of  roll  angle  upon  the  natural  convection  heat 
transfer  from  an  array  of  cylindrical  spines  or  pin  fins.  A 
correlation  is  presented  that  permits  evaluation  of  the 
Nusselt  Modulus  as  a  function  of  the  Rayleigh  Modulus  and 
the  roll  angle  and  it  shows  that  maximum  heat  transfer  is 


obtained  at  roll  angles  between  20  and  30  degrees.  Kt 
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I.  INTRODUCTION 


z'  Electronically  steered  phased  array  antennae  frequently 
exhibit  degraded  performance  due  to  the  lack  of  proper  heat 
transfer  between  the  interior  crystalline  materials  and  the 
environment.  This  causes  higher  than  optimum  temperature 
environments  for  the  crystalline  materials  and  a  subsequent 
loss  of  reliability^-*-}^ 

One  possible  solution  to  such  a  problem  considers  a 
surface  augmentation  through  the  use  of  pin  fins  (or  spines) 
to  assist  in  the  dissipation  of  the  heat.  -'-/■{!<■  1  *'  ' 

The  problem  is  more  severe  when  the  installation  is 
stationary  so  that  the  dissipation  from  the  spines  is  in  the 
natural  convection  mode  and,  because  of  this,  the  focus  of 
this  investigation  is  in  the  area  of  natural  convectio*.  heat 
transfer  and  methods  of  enhancement  in  this  domain. 

Convective  heat  transfer  from  horizontal  cylinders  has 
been  the  subject  of  extensive  studies  which  were  summarized 
by  Morgan[2].  In  addition,  studies  of  heat  exchanger  tube 
bundles  are  summarized  by  Incropera  and  and  DeWitt[3]. 
However,  after  a  considerable  search  of  the  literature,  only 
a  small  amount  of  experimental  data  and  analytical  correla¬ 
tions  dealing  with  natural  convection  in  an  array  of  hori¬ 
zontal  cylinders,  has  been  found. 

Intuitively,  one  feels  that  the  orientation  of  the  array 
must  have  an  effect  on  the  heat  transfer  but  apparently  no 
known  experimental  correlations  have  been  formulated.  Thus, 
the  focus  of  this  thesis  is  in  the  effect  of  roll  angle  on 
the  heat  transfer  rates  from  an  array  of  cylindrical  spines. 

Additional  impetus  for  investigations  in  similar  areas 
with  cylindrical  spine  arrays  derives  from  the  electronics 
industry  where  such  cooling  arrangements  are  used  in  the 
external  cooling  of  dense  integrated  circuit  arrangements. 


Dimensionless  correlations  using  dimensionless  parame¬ 
ters  are  usually  developed  for  presentation  of  experimental 
results  and  this  allows  for  possible  application  to  other 
systems.  In  developing  the  solution  for  the  heat  transfer 
across  a  boundary  layer,  it  has  been  shown  by 
Schlichting[4 ],  and  others  [3,5],  that  the  Nusselt  number 
(Nu)  is  a  proper  dimensionless  parameter  for  natural  or 
forced  convection.  In  forced  convection,  the  Nusselt  number 
is  a  function  of  the  Reynolds  number  (Re)  and  the  Prandtl 
number  (Pr)  .  However,  in  pure  natural  convection  where  the 
buoyancy  is  the  only  driving  fcrce,  the  fluid  velocity  is 
determined  entirely  by  the  quantities  contained  in  the 
Grashof  number  (Gr)  which  is  used  instead  of  the  Reynolds 
number  to  establish  dynamic  similarity.  The  Grashof  number 
is  developed  quite  thoroughly  by  Kreith  and  others  [3,5]  and 
guickly  appears  in  a  dimensional  analysis. 

In  addition,  for  gases  the  Prandtl  number  (Pr)  is  nearly 
constant  so  that  the  Nusselt  Modulus  (Niy)  can  be  considered 
to  be  a  function  of  the  product  of  the  Grashof  number  and 
the  Prandtl  number.  This  product  is  defined  as  the  Rayleigh 
or  number  (Ra)  and  is  discussed  to  further  detail  in  the 
following  paragraphs. 

In  order  to  establish  the  heat  transfer  coefficient  and 
to  develop  a  correlation  of  the  data  with  dimensionless 
parameters,  relationships  for  the  heat  transfer  coefficient 
and  the  Nusselt  number  were  established  and  measurements 
were  taken  which  permitted  its  calculation.  A  summary  of 
the  derivation  is  provided  in  Figure  1.1  and  a  detailed 
rivation  is  provided  in  Appendix  A.  The  relation  holds 
for  the  condition  that  no  heat  transfer  occurs  at  the  fin 
tip.  Such  a  condition  can  be  effectively  achieved  by  insu¬ 
lating  the  tip. 

After  calculating  the  heat  transfer  coefficient  the 
Nusselt  number  was  then  determined  from  the  relation: 


APPENDIX  A  > 

DEFINITION  OF  THE  HEAT  TRANSFER  COEFFICIENT 

Fins  and  spines  (pin  fins)  of  various  geometries  and 
thermal  conductivities  respond  differently  to  identical  and 
uniform  heat  sources  and  sinks.  Similarly,  there  are 
numerous  ways  in  which  the  temperatures  and  heat  transfer 
coefficients  of  sources  and  sinks  may  vary.  Important  to 
the  analysis  of  fin  geometries  are  the  constraints  or 
assumptions  which  are  employed  to  define  and  limit  the 
problem  and  often  to  simplify  its  solution.  The  analysis  of 
fins  and  spines  employ  the  assumptions  of  Murray[7]  and 
Gardner[8].  These  limiting  assumptions  are: 

1.  The  heat  flow  in  the  fin  and  its  temperatures  remain 
constant  with  time. 

2.  The  fin  material  is  homogeneous,  its  thermal  conduc¬ 
tivity  is  the  same  in  all  directions  and  it  remains 
constant. 


3.  The  heat  transfer  coefficient  between  the  fin  and 
the  surrounding  medium  is  constant  and  uniform  over 
the  entire  surface  of  the  fin  (spine) . 

4.  The  temperature  of  the  medium  surrounding  the  fin  is 
uniform. 


5.  The  fin  thickness  is  so  snail  compared  with  its  height 
that  temperature  gradients  across  the  fin  thickness 
may  be  neglected. 

6.  The  temperature  at  the  base  of  the  fin  is  uniform. 


7.  There  is  no  contact  resistance  where  the  base  of  the 
fin  joins  the  prime  surface. 


V.  CONCLUSION  AND  RECOMMENDATIONS 


A.  CONCLUSION 

The  results  of  the  experimental  program  indicated  that  a 
dependence  on  roll  angle  did  exist.  This  dependence  is 
relatively  minor  (i.e.,  on  *-he  order  of  3%),  however  for 
temperatures  in  the  range  of  200  degrees  Fahrenheit  this 
could  mean  a  difference  of  5  to  6  degrees  in  the  surface 
temperature  of  an  operating  system.  Therefore,  further 
investigation  at  higher  Rayleigh  numbers  may  be  deemed 
appropriate. 

B.  RECOMHENDATIOHS 

The  following  recommendations  are  provided  for  possible 
follow  on  projects  of  a  similar  nature. 

1.  Conduct  the  experimental  pregram  with  the  goal  of  deter¬ 
mining  and  removing  the  radiation  heat  transfer. 

2.  Investigate  other  devices  and  methods  of  temperature 
measurement  in  order  to  remove  more  of  the  uncertainty 
associated  with  the  use  of  thermocouples. 

3.  Conduct  a  similar  experiment  for  the  pitch  angle. 

4.  Attempt  a  similar  program  with  water  as  the  convecting 
medium. 


which  correlated  experimental  results  of  several  researchers 
very  well. 

The  values  obtained  in  this  experimental  program  agreed 
favorably  with  these  predictions  with  an  average  of  19% 
deviation  throughtout  the  range  investigated.  Noting  that 
the  effect  of  heat  transfer  by  radiation  were  also  present, 
the  deviation  is  attributed  to  radiation  losses  primarily. 
These  losses  would  appear  to  net  have  a  detrimental  effect 
on  the  observed  roll  angle  dependence. 


The  values  or  Z  ($)  are  tabulated  ia  Table 


1. 


TABLE  1 


Enhancement 

Factor 

Roll 

Enhancement 

Angle(0) 

Factor  (Z(0)) 

0° 

0.5916 

10° 

0.5986 

20° 

0.6075 

30° 

0.5990 

45° 

0.5928 

The  values  of  Nusselt  numbers  obtained  from  the  experi¬ 
mental  results  are  a  function  of  a  given  Rayleigh  number  and 
roll  angle  and  may  be  correlated  via: 


Nud 


Z(0)Ra 


0.233 


(4.1) 


However,  observe  that  the  dependence  on  the  roll  angle  is 
weak  (Table  1)  and  no  mathematical  relationship  was  deter¬ 
mined  for  the  values  of  Z  (#)  . 

Results  were  also  compared  with  several  previously 
published  correlations  for  the  single  horizontal  cylinder. 
Churchill  and  Chu[t>]  suggested  an  equation  of  the  fora 


f 


Nu 


*  V 


.60  + 


0.387  Ra 


(1/6) 


1  + 


[0.559 
,  Pr 


,(9/16) 


(8/27) 


IV.  BE5DLTS  AMD  COE RELATION 


A.  RESULTS 

During  each  data  collection  run  the  temperatures  at  the 
base  and  tip  of  each  spine,  as  veil  as  four  environmental 
temperatures  were  recorded.  These  temperatures  were  used  as 
inputs  for  calculation  of  the  various  output  parameters  for 
each  spine.  These  sixteen  calculations  were  then  averaged 
to  provide  the  data  point  for  that  set  of  conditions.  Forty 
seven  data  points  were  obtained.  Typical  data  runs  are 
included  as  Appendix  C. 

At  a  specified  angle,  the  natural  logarithm  of  the 
Nusselt  number  was  plotted  as  a  function  of  the  natural 
logarithm  of  the  Rayleigh  number  for  each  of  the  data  points 
at  a  specified  angle.  These  plots  are  included  as  Figures 

B.  1  through  B.5  in  Appendix  B.  This  graphical  display  was 
then  used  in  order  to  facilitate  correlation  of  the  test 
data . 

The  plotted  results  indicated  that  a  small  dependence  on 
roll  angle  did  exist.  Figures  B.6  and  B.7  were  plotted  to 
show  these  effects  more  specifically.  Figure  B.7  was 
plotted  only  for  Rayleigh  numbers  between  579  and  586 
because  of  the  narrowness  of  the  overall  Rayleigh  number 
range  obtained. 

B.  CORRELATION 

For  each  of  the  curves  showing  log  {Nu)  as  a  function  of 
log  (Ra) ,  a  least  sguares  curve  fit  routine  was  used  and  a 
linear  relation  established  tc  approximate  these  curves 
following  a  method  presented  by  Kraus  and  Bar-Cohen[ 1  ].  A 
correlating  equation  was  obtained  for  each  of  the  roll 
angles  (0) ,  where  Z  (0)  defines  a  function  of  the  roll  angle. 


B.  DATA  BEDUCTIOH 


Temperature  data  printed  by  the  Autodata  Nine  during  the 
individual  runs  was  used  as  input  parameters  (corrected  from 
the  calicration  curve  for  each  thermocouple)  for  calculation 
of  the  desired  parameters  for  plotting  and  evaluation. 

All  desired  output  parameters  were  calculated  by  the 
computer  using  appropriate  relationships  for  the  thermophy¬ 
sical  properties  £3]  as  functions  of  temperature  and  meas¬ 
ured  temperatures  as  input  paraneters.  The  computer  program 
used  for  the  calculations  is  provided  in  Appendix  E. 

Graphical  display  of  the  data  was  also  utilized  to  aid 
in  data  analyses.  These  curves  are  shown  as  Figures  B.  1 
through  B.7. 


III.  EXPERIMENTAL  PBOCEDPBE 


A.  DATA  GATHEBI1G 

The  experimental  procedure  was  to  set  the  disk  at  the 
desired  roll  angle  using  a  liguid  bubble  level  indicator  and 
protractor  to  measure  the  angle.  The  temperature  control 
thermostat  of  the  environmental  chamber  was  then  set  to  the 
desired  temperature.  Steady  state  environmental  conditions 
were  attained  in  approximately  six  to  eight  hours.  Next  the 
heaters  imbedded  in  the  disk  were  energized  to  elevate  the 
temperature  of  the  disk  and  base  of  the  spine  to  a  tempera¬ 
ture  chosen  for  a  predetermined  Rayleigh  number.  The  input 
voltage  for  the  heaters  was  varied  and  the  temperatures 
monitored  until  a  steady  state  condition  was  reached.  This 
was  indicated  by  constant  input  power  and  constant  tempera¬ 
tures.  The  time  to  attain  a  total  steady-state  condition 
was  approximately  three  to  five  hours  in  addition  to  the 
time  reguired  for  the  chamber  tc  stabilize. 

A  steady  state  condition  was  assumed  to  exist  when  no 
change  in  temperature  occurred  for  four  consecutive  tempera¬ 
ture  readings  taken  at  fifteen  minute  intervals.  These 
temperatures  were  then  used  for  the  input  data  for  the 
calculations.  Typical  dataruns  are  presented  in  Appendix  B. 

After  completion  of  data  taking  at  a  set  condition, 
changes  were  made  in  the  roll  angle,  the  temperature  of  the 
system  or  temperature  of  the  environment,  and  the  process  of 
attaining  and  maintaining  a  new  condition  of  steady  state 
aain  started. 


was  a 


figure  2.4  lest  Disk  and  Spine  Arrangenent 


E.  ENVIRONMENTAL  CHAMBER 

The  test  system  was  housed  in  an  environmental  chamber 
in  which  the  ambient  temperatures  could  be  maintained  over  a 
range  between  35  and  135  degrees  Fahrenheit. 


Figure  2.3  Sketch  of  Front  View  of  Test  Disk 
D.  ACDREX  AUTODATA  NINE 

The  thermocouple  modules  in  the  data  acquisition 
system’s  thermocouple  modules  include  an  ice-point  compensa¬ 
tion  network  and  provide  automatic  data  scanning  and 
printout  for  up  to  100  channels  of  input.  Due  to  the 
internal  ice-point  compensation  networks,  no  external  ice 
baths  or  reference  junctions  were  required.  In-place  cali¬ 
bration  checks  were  made  at  various  temperatures  over  the 
expected  temperature  range  in  order  to  develop  temperature 
correction  factors  for  each  thermocouple,  its  connecting 
cabling  and  the  Autodata  Nine  circuitry  as  a  unit,  in  addi¬ 
tion,  the  complete  disk  and  spine  array  was  immersed  in  an 
ice-bath  in  order  to  establish  one  of  the  check  points. 


Figure  2.2  Test  Disk  and  Spine  Assembly 


spine)  to  the  outer  level  (exterior  surface)  of  the 
insulation . 

C.  THERHOCOOPLES 

All  thermocouples  were  fabricated  from  copper-constantin 
(Type-T)  of  AWG  No.  30  wire.  They  were  butt  welded  with  a 
spherical  bead  of  approximately  0.05  cm  diameter  and 
smaller.  All  interconnecting  cables  were  also  of  copper- 
constantin  in  order  to  minimize  the  uncertainty  and  to 

preclude,  as  much  as  possible,  spurious  thermoelectric 

effects.  All  thermocouples  were  connected  to  the  ACUREX 

Autodata  Nine  data  acquisition  system  for  temperature 

recording. 


II.  EXPERIMENTAL  APPARATUS 


Figure  2. 1  displays  the  configuration  of  the  experi¬ 
mental  apparatus  and  the  techniques  used  in  data  acquisi¬ 
tion.  The  following  paragraphs  provide  a  more  indepth 
description  of  the  component  parts. 

A.  TEST  DISK 

The  disk,  shown  in  Figures  2.2  and  2.3,  was  fabricated 
entirely  of  2024-T4  aluminum  plate, 6.35  cm  in  thickness.  A 
circular  disk  of  diameter  equal  to  22.86  cm  was  cut  from 
this  plate  and  fitted  with  an  array  of  16  cylindrical  spines 
as  the  heat  transfer  surfaces.  Six  equally  spaced  radial 
holes  were  drilled  around  the  perimeter  of  the  disk  for 
insertion  of  six  200  watt  cartridge  heaters  which  provided 
the  heat  input  necessary  to  elevate  the  temperature  of  the 
spines.  The  heaters  were  powered  by  a  controllable  power 
source  in  order  to  attain  a  steady-state  condition  for  data 
measurements.  The  disk  was  fully  insulated  with  soft  rubber 
refrigeration  insulation  (kx  *  0.  15  w/m-K)  to  minimize  heat 
losses  from  the  disk  surface  due  to  both  radiation  and 
convection. 

B.  CYLINDRICAL  SPINES 

The  spines,  shown  in  Figures  2.3  and  2.4  were  similarly 
constructed  of  the  2024-T4  aluminum.  Two  sizes  of  spines 
were  used  in  data  gathering  (0.9525  cm  diameter  x  3.75  cm 
length  and  1.27  cm  diameter  x  5.08  cm  length)  to  provide  a 
wider  range  of  Rayleigh  Numbers.  Extension  pieces,  also 
made  of  the  same  aluminum,  were  necessary  to  extend  the 
actual  base  of  the  heat  transfer  surface  (the  cylindrical 


Ra  =  GrPr 


(1.2) 


where  again  the  subscript  "d"  infers  a  dependence  on  the 
spine  diameter  and  where  the  symbols  are  as  defined  in  the 
nomenclature  list. 


1  3 


Figure  1.1  A  Gelation  for  the  Heat  Transfer  Coefficient 


(1.1) 


where  Nu^  is  the  Nusselt  number  based  on  the  spine  diameter. 

The  Bayleigh  number  was  determined  from  the  product  of 
the  Trash of  number  and  the  Prandtl  number  which  were  calcu¬ 
lated  from  the  following: 


Gr 


g  B(T  -  T  )  d“ 
s 


C  U 

Pr  =  -T2” 
Kf 
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8. 


There  are  no  heat  sources  within  the  fin  itself 


9.  The  heat  transferred  through  the  outermost  edge  of 
the  fin  is  negligible  compared  with  that  leaving  the 
fin  through  its  lateral  surface. 

10.  Heat  transfer  to  or  from  the  fin  is  proportional  to 
the  temperature  excess  between  the  fin  and  the  sur¬ 
rounding  medium. 

In  addition  to  proposing  a  profile  function  for  longitu¬ 
dinal  fins  Gardner£8]  proposed  a  profile  function  for 
spines: 


f2(x)  = 


(A.1) 


With  the  appropriate  value  of  n,  eguation  A.1  may  be  used 
for  the  development  of  a  generalized  differential  eguation 
for  spines. 

Figure  A.1  shows  a  spine  of  arbitrary  profile.  It  can 
be  seen  that  the  spine  cross  section  normal  to  the  flow  of 
heat,  the  confining  profile,  and  the  perimeter  of  the  spine 
are  perfectly  arbitrary  functions  of  the  distance  x  from  the 
tip  of  the  spine.  A  differential  eguation  for  the  tempera¬ 
ture  excess,  0  =  t  -  t  ,  may  he  written  by  considering  the 
heat  conduction  into  and  out  of  the  element  dx  through  the 
cross  section  f  (x)  : 


A 

--r—f 

\  i 


,P=£(x) 


-jr6M 


Figure  A. 1  Spine  of  Arbitrary  Profile 

The  difference  in  heat  flow  into  and  out  of  the  element 
must  equal  that  dissipated  by  the  surface  of  the  spine.  If 
dissipation  occurs  by  convection  and  h  is  the  convection 
coefficient 

dq  =  h  fg(x)  0  dx 


where  f3(x)  defines  a  perimeter  function  P  ( x)  which  depends 
upon  the  distance  x  from  thn  origin  of  the  coordinate 
system.  When  convection  and  conduction  are  equated, 


4s[f2(!,>3!]= 


hf  3 ( x)0 


Upon  rearrangement  the  qonecalized  differential  equation 
becomes 


<re  df,  (X)  dQ 


ij 


% ; 


The  relationship  between  f,  (x)  and  £2(x)  is 


rx(x)  =  Ti  |f2(x)J' 


where  f2(x)  was  defined  by  equation  A.  1.  Thus  equation  A. 2 
may  be  written 


[f2(xf 


dx  dx 


d  ff,(x)l2d0 

+  —  - 

/4..  *  J.. 


2h  f2Cx)0  (A. 3) 


dx  k 


where 


f3^x)  =  2nf?(x) 


Equation  A.  3  is  a  second-order  differential  equation 
with  variable  coefficients  except  when  the  spine  cross 
section  normal  to  heat  flow  is  constant.  It  nay  be  solved 
by  termwise  comparison  with  the  generalized  Bessel  equation. 

For  the  cylindrical  spine,  shown  in  Figure  A. 2 


f2(x) 


3s[f2(*>]: 


For  this  case  the  exponent  in  equation  A.1  has  a  value  of 
zero  correspond ii.  j  to  n  =  1/2 .  Actually  £  in  equation  A.1 

is  replaced  by  d,  the  spine  diameter.  Using  this  value  of  f- 
the  general  di f f erent ia 1  equation  for  the  temperature 


excess  tecoiaes 


Figure  A. 2  Cylindrical  Spine 


2 

m0  = 


0 


where  m  =  (ah/kd?  . 

The  general  solution  to  equation  A. 4  is 


9(x)  =  C,emx  +  c  e-mx 
1  ^2 


where  the  arbitrary  constants  are  determined  from 
boundary  conditions 


This  then  leads  directly  to  the  equation  for  the  temperature 
excess  profile 


e(x)  =  qk  ggsjijg* 
b  cosh  mb 


(A. 7) 


and  the  heat  flow  entering  the  fin  (spine)  base  is 


tt  ,  ,2  Q 
— ; — kd  m6, 
4  b 


tanh(mb) 


(A. 8) 


The  value  of  h  may  be  determined  from  equation  (A. 7)  by 
taking  temperature  readings  at  the  tip  and  the  base  which 
are  precisely  "b"  meters  apart. 

Thus  with  base  and  tip  temperature  excesses, 

«(x=b)  ebe-(x=a)  =  ea 
It  is  easy  to  see  from  equation  (A. 7)  that 

e  .  _!k_ 
a  cosh  mb 


or  that 


cosh ( mb ) 


from  which 


mb  =  arccoshCQ,  /9  )  =  a 
b  a 


(A.  9) 


Then  with  a 


(4h/kd 


one  may  then  solve  for  h. 


kd,a,.  2 


Fijure  B.2  logfiNu  vs.  lojeRa  (10  degrees  roll). 
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ROLL  ANGLE  (0)  (DEGREES) 


APPENDIX  C 


TYPICAL  DATA  RUNS 
TABLE  C.l 


Typical  Results  (579  <  Ra  <  586) 


.89  0  = 

0° 

U 

T.  ( °  K  ) 

D 

T  (  °  K ) 
o 

Ra 

Nu 

Nu  , 
ch 

n 

( w/m^  ‘  ° 

312.48 

312.25 

579  .  40 

2.2339 

2.3333 

6  .162 

312.32 

312.05 

579.31 

2 .5  850 

2.3283 

7.128 

312.40 

312.23 

579.40 

2.9536 

2.3333 

7.043 

312 . 32 

312.03 

572.94 

2.8039 

2.3280 

7.731 

312.37 

312.08 

5  74  .  69 

2.7417 

2.3295 

7.560 

312.59 

312  .33 

582  .4  3 

2.4842 

2.3358 

6.853 

312  .61 

312.35 

583.17 

2.4806 

2.3364 

6.843 

312.64 

312.34 

583.84 

2 . 7984 

2.3367 

7.720 

312.41 

312.16 

576.73 

2 .4066 

2  .3311 

6.637 

312.44 

312.16 

577.28 

2.7274 

2.3316 

7.522 

312.58 

312.31 

581.97 

2.6458 

2.3354 

7.298 

312.66 

312  .38 

584.27 

2.5824 

2.3373 

7.124 

312.61) 

312.3? 

582.43 

2.6995 

2.3358 

7.446 

312.  3 7 

312.14 

575.71 

2.2505 

2.330  3 

6  .  206 

n  2 . 4  4 

312.26 

680.31 

2.7625 

2.3340 

7.620 

3  1  2  .  4  8 

312.19 

578.48 

2.7727 

2.3326 

7.647 

.9  3  0  r 

It)0 

312. 76 

312.47 

586.28 

2 .7304 

2.3389 

7.534 

312.69 

312.33 

581.14 

2.5429 

2.3347 

7.015 

312.73 

312.51 

587.20 

2.5123 

2.3396 

6.932 

3  1  2  .  5  9 

312.31 

580.77 

2.7581 

2.3344 

7.609 

312.64 

312.36 

582.51 

2.6972 

2.3358 

7.441 

312.81 

312.55 

588.37 

2.4534 

2  .  3406 

6  ,  :  7  0 

312.89 

312.63 

590  .9  5 

2.4434 

2.3426 

6.743 

312.86 

312.57 

589.49 

2.7638 

2.3415 

7.627 

312.63 

312.38 

582.70 

2.3763 

2.3360 

6.556 

312.67 

312.38 

583.25 

2.6932 

2.3364 

7.430 

312.86 

312.58 

589.76 

2.6057 

2.3417 

7.191 

312.87 

312.60 

568.37 

2.5903 

2.3420 

7.038 

312.82 

312.54 

590.21 

2.6658 

2 .3406 

7  .356 

312.59 

312.31 

580.77 

2.7581 

2.3344 

7.609 

312.76 

312.47 

580.25 

2.7304 

2.3389 

7.534 

312.71 

312.42 

5  84  .  4  4 

2.7380 

2.3  374 

7.554 

» 

TABLE  C.l 

i 

Typical 

Results 

(579  <  Ra  <  586) 

( cont ' d . ) 

1 

s 

=  296.93  0 

=  20° 

■-V  i 

312 . 71 

312.42 

584.44 

2 .7380 

2 .3374 

7.554 

2 

312.54 

312.27 

579.29 

2.5523 

2  .  3332 

7.040 

3 

312 . 72 

312.46 

585.36 

2 .5216 

2.3381 

6.958 

Q  '* 

312 . 54 

312.25 

578.92 

2.7684 

2.3329 

7.636 

312  .  59 

312 . 31 

580.68 

2.7047 

2.3343 

7.461 

312.76 

312.44 

585 . 84 

2.9965 

2 .3384 

8.268 

7 

312.78 

312 . 52 

587.28 

2.4811 

2.3397 

6.791 

8 

312.81 

312.51 

587.65 

2.7738 

2  .  3400 

7.654 

0 

312.58 

312  .33 

580 .86 

2 . 3851 

2  .  3345 

6.579 

r  ^ 

312.56 

312.27 

579.57 

2 . 7130 

2.3343 

7.484 

i  i 

312.75 

312.47 

586.09 

2 .6246 

2 .3387 

7.242 

1 2 

312.82 

312 . 54 

588.37 

2.5597 

2 . 3406 

7.063 

i  3 

312.77 

312.48 

586.56 

2.6733 

2.3391 

7.376 

1  '1 

312 . 54 

312.25 

578.92 

2.7680 

2  .  3329 

7.636 

;'v;  i  s 

312.59 

312.31 

580.77 

2.7584 

2.3344 

7.609 

a  ir- 

312.69 

312  .  31 

580.77 

2.7584 

2.3344 

7.609 

*  •  T* 

>v  S 

=  296.94  0 

=  30° 

::v  i 

312.65 

312 . 36 

582 . 81 

2 . 7480 

2 .3359 

7.581 

O*  ? 

312.54 

312.27 

579.29 

2.5523 

2.3332 

7.040 

a 

312.67 

312.40 

583.52 

2.5332 

2 .3366 

6.989 

* 

312  .  64 

312.25 

578.92 

2 . 7684 

2.3329 

7.636 

312.53 

312.31 

579 . 76 

2.1768 

2.3336 

6.005 

», 

312. 7C 

312.44 

5  8  5.64 

2.9965 

2  .  3  384 

8.268 

:  v  7 

312. 7 6 

312.52 

587.28 

2.4611 

2.3  397 

6.791 

h 

312.58 

312.27 

679 .94 

2.9250 

2.3337 

8.068 

312.56 

312.27 

579  .97 

2.7130 

2  .3334 

7.484 

10 

312.70 

312.51 

586.73 

2 .2509 

2.3392 

6.211 

1  i 

312  .69 

3^1  2.42 

584  .27 

2 .6842 

2.3373 

7.268 

1 2 

312.62 

312.24 

588.37 

2.6599 

2  .  3406 

7.063 

1  .5 

312.77 

312.48 

586.56 

2.6733 

2 .3391 

7.376 

■ ;/  ]  4 

312.48 

312 . 14 

575.25 

2  .  7890 

2.3299 

7.692 

pa  ^  ^ 

312.59 

312.31 

580.77 

2.7581 

2 .3344 

7.609 

1 

312  .  59 

312.31 

580.77 

2.7581 

2.3344 

7.609 

1 

40 

• 

| 

.■■7  ■  2  /’A 

afcf .  *  -  a  .  -  * 

»  •  ■,  *,  •  ,  **'  1 
•'./•/•A'-' 

TABLE  C.l 


Typical  Results  (579  <  Ra  <  586)  (cont'd.) 


296.96  0  =  45° 


312.65 

312.42 

582.12 

312.48 

312.22 

576.06 

312.67 

312 .40 

582.12 

312 . 48 

312 . 19 

575.69 

312.53 

312.25 

577.44 

312 . 76 

312.49 

585.15 

312 .78 

312 . 52 

585 . 89 

312.75 

312 .51 

585.35 

312.58 

312.33 

579.47 

312.50 

312 .22 

576.34 

312.69 

312.42 

582.86 

312.82 

312.54 

586.99 

312.71 

312.43 

583.31 

312.43 

312.14 

573.85 

312.59 

312.31 

579.38 

312.59 

312.31 

579.38 

2.21640 

2 .33550 

6.115 

2.56722 

2 .33057 

7.081 

2 .53847 

2.33550 

7.004 

2 .78452 

2 . 33024 

7.681 

2.72293 

2 .33169 

7.511 

2.46985 

2 .33796 

6 .815 

2.46626 

2 .33856 

6  .  806 

2.25553 

2 .33812 

6.224 

2 . 39009 

2.33335 

6.594 

2.72653 

2 . 33079 

7.521 

2.63975 

2 .33611 

7.283 

2.56502 

2 .33945 

7.076 

2.69100 

2 .33647 

7.425 

2 . 79496 

2 . 32875 

7.709 

2 . 76402 

2.33328 

7.625 

2.76402 

2 . 33328 

7.625 

TABLE  C.2 


Typical  Results  (1305  <  Ra  <  1380) 
317.83  0  =  0° 


h 


Tb(°K) 

Tq(°K) 

Ra 

Nu 

Nuch 

(w/m^ ‘ 

338.49 

338.04 

1304.18 

2.97416 

2.76085 

6.553 

338.49 

337.99 

1302.95 

3.27379 

2.76029 

7.213 

338.52 

338.06 

1305.56 

3.07919 

2.76146 

6.784 

338.44 

337.04 

1300.20 

3 .28048 

2.75905 

7.227 

338.53 

338.05 

1305.56 

3.15390 

2.76146 

6.949 

338.49 

338.02 

1303.57 

3.12388 

2.76057 

6.882 

338.52 

338.05 

1305.40 

3.11699 

2.76139 

6.868 

338.50 

338.01 

1303.72 

3.23238 

2.76064 

7.122 

338.53 

338.08 

1306.47 

3.00352 

2.76188 

6.618 

338.53 

338.08 

1306 . 31 

2.96674 

2.76180 

6.537 

338.07 

338.21 

1313.81 

2.98347 

2.76517 

6.575 

338.03 

338.17 

1311.66 

3.06200 

2.76421 

6.747 

338.50 

333.01 

1303.73 

3.23238 

2.76064 

7.122 

338.51 

338.04 

1304.78 

3.04564 

2.76111 

6.710 

338.55 

338.07 

1306.63 

3.18820 

2.76195 

7.025 

338.51 

338.01 

1304.02 

3 . 30642 

2.76077 

7.285 

314.72  0  = 

10° 

334 . 82 

334.35 

1326.67 

3.23972 

2 .77102 

7.075 

334 . 87 

334  .  42 

1330.20 

3.07632 

2.77259 

6.718 

334.81 

334.33 

1325 . 54 

3.28189 

2.77052 

7.166 

334.79 

334 .35 

1325.88 

3.04951 

2.77067 

6.659 

334.76 

334.32 

1323.93 

3.97872 

2 . 76980 

6 . 504 

334 . 82 

334.38 

1327.63 

3.00764 

2.77145 

6,568 

334 .75 

334.28 

1322.49 

3.24195 

2.76916 

7.101 

334 .88 

334.39 

1329.72 

3 . 34543 

2.77238 

7 . 306 

334 .87 

334.42 

1330.06 

3.03809 

2.77253 

6.635 

334.81 

334.37 

1326 . 84 

3.05689 

2.77110 

6.654 

334 . 88 

334.41 

1329.89 

3 .23042 

2 .77245 

7.055 

3  34.92 

334.41 

1331.16 

3.45717 

2.77302 

7.550 

334.95 

3  34  .  4  9 

1334.06 

3.06572 

2.77430 

6.696 

334 . 80 

334.34 

1325.88 

.12721 

2.77067 

6.829 

334.93 

334.47 

1333.24 

3.18207 

2 .77394 

6.950 

334 . 80 

334.33 

13  25.54 

3.20581 

2 .77052 

7.000 

TABLE  C . 2 


Typical  Results  (1305  <  Ra  <  1380)  ( Cont . ) 


=311.24  0  =  10° 


330 . 59 

330.36 

1347.11 

3.11643 

2.78019 

6.736 

330.48 

330.07 

1340  .  81 

3.01289 

2 .77741 

6.512 

330.48 

330.04 

1339.97 

3.21805 

2 .77704 

6.955 

330.43 

329.98 

1336.56 

3 . 22605 

2.77553 

6.972 

330.59 

330.14 

1346.60 

3.23823 

2 .77997 

7.000 

330.54 

330.07 

1342.69 

3 .37037 

2 .77824 

7.285 

330.43 

330.00 

1337.07 

3.10318 

2.77576 

6 . 706 

330.59 

330.16 

1347.11 

3.11653 

2  .  78019 

6.736 

330.54  • 

330.09 

1343.38 

3.20669 

2 .77855 

6.931 

330.43 

329.98 

1336.56 

3.22605 

2  .  77553 

6.972 

3  30 . 59 

330.14 

1346.60 

3.23823 

2.77997 

7.000 

330.60 

330 .14 

1346.60 

3.31969 

2 .77997 

7.176 

330.53 

330.08 

1342.87 

3.25068 

2.77832 

7.026 

330.48 

330.05 

1340 . 30 

3.13528 

2.77719 

6.776 

330 . 59 

330 . 16 

1346 .93 

3.15756 

2.78012 

6 . 825 

330 .37 

329.92 

1333.00 

3.27927 

2.77396 

7.086 

310.65  0  = 

20° 

330 . 32 

329.88 

1377.21 

3.06971 

2.79335 

6.628 

330.32 

329 . 89 

1377.54 

2.98876 

2.79350 

6.453 

330.26 

329.82 

1373.47 

3.16018 

2.79173 

6.823 

330 . 26 

329 . 82 

1373.47 

3.16018 

2.79173 

6.823 

330 . 32 

329 . 87 

1376.70 

3.18968 

2.79313 

6 .887 

330.26 

329 . 80 

1372 .95 

3.28058 

2.79151 

7.082 

330.26 

329 . 78 

1372.26 

3 .44194 

2 .79121 

7.431 

330.37 

329 .88 

1378.91 

3.46225 

2 . 79409 

7.476 

330.43 

329.93 

1381.96 

3.53134 

2.79541 

7.625 

330 . 32 

329.87 

1376 .88 

3.14908 

2.79321 

6.799 

330.43 

329.98 

1383.67 

3.13060 

2.79614 

6 . 760 

330.43 
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330.43 

329.98 
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2.79614 
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330.26 

329 . 82 
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9 
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10 
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TABLE  C .  2 


Typical 

Results 

(  1305  <r 

Ra  <  1380)  (Cont.) 

-  .  4 

hr--: 

T  = 

c 

301 

.86  0  = 

30.0° 

vZv 

o 

1 

318.91 

318.58 

1378.97 

2 .79324 

2 . 79445 

5.869 

Ur? 

2 

318 . 87 

318.48 

1373.41 

3.27655 

2 . 79204 

6.884 

•V-V 

3 

318.96 

318.56 

1379 .78 

3.35353 

2 . 79480 

7.047 

4 

318.96 

318.57 

1379.99 

3.21071 

2.79489 

6.747 

r 

.j 

319.00 

318.62 

1383.17 

3.20213 

2.79627 

6.729 

*v*\: 

6 

318.97 

318 . 59 

1381.18 

3.20759 

2 . 79540 

6.740 

r 

7 

319.02 

318.63 

1384.15 

3.24665 

2.79669 

6.823 

K  -  v 

8 

318.95 

318.56 

1379.36 

3.25983 

2.79462 

6.850 

' .  **» 

9 

318.89 

318.52 

1375.58 

3.12968 

2 .79298 

6.576 

/•y« 

10 

318.94 

318.57 

1379.57 

3.11910 

2.79471 

6  .554 

2] 

318.95 

318.57 

1379.5? 

3.21197 

2.79471 

6.749 

,*•  .**< 

]  2 

319  .01 

318.58 

1382.16 

3 .53414 

2.79583 

7.426 

f 

13 

318.99 

318.56 

1380 . 76 

3.58582 

2.79522 

7.535 

14 

318.96 

318 . 59 

1380.76 

3.02125 

2.79522 

6  .  349 

15 

318.91 

318.54 

1377.40 

3.07746 

2 . 79377 

6.466 

•vV 

*  V*  - 

16 

318.98 

318.59 

1381.35 

3.25438 

2.79548 

6.838 

•T.\' 

c 

V 

309 

.50  0  = 

4  5° 

1 

328.21 

327.79 

1339 . 79 

3.07511 

2.77703 

6.6]] 

2 

328 . 26 

3  2  7.83 

1342.74 

3.19372 

2.77833 

6 . 866 

3 

328.27 

327.87 

1343.98 

2.97942 

2 . 77887 

6.406 

4 
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APPENDIX  D 


DATA  REDUCTION  PROGRAM 
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*  PROGRAM  TO  CALCULATE  HEAT  TRANSFER  COEFFICIENT  ON  A  * 

*  PINNED  CIRCULAR  DISK  * 

****#**###*###**##*##*#*###*#*##**##*##***#**####*#***#*## 

##***#****#*****#DgFjf! j ij> jon  OF  VARIABLES****************** 
A  =  FACTOR  IN  CHURCHILL  AND  CHU  CORRELATION  USED  HERE 
TO  SIMPLIFY  AND  SHORTEN  THE  EXPRESSION 
B  =  LENGTH  OF  THE  SPINE 

C  =  SIMPLIFYING  FACTOR  FOR  CHURCHILL  AND  CHU  CORRELA¬ 
TION 

CB  =  CORRECTION  TO  THE  SPINE  BASE  TEMPERATURE 

CO  =  CORRECTION  TO  THE  SPINE  TIP  TEMPERATURE 

D  =  DIAMETER  OF  SPINE 

E  =  SIMPLIFYING  FACTOR  FOR  THE  CHURCHILL  AND  CHU  COR¬ 
RELATION 

G  =  GRAVITATIONAL  ACCELERATION 

GR  =  GRASHOF  NUMBER 

H  =  HEAT  TRANSFER  COEFFICIENT  (EXPERIMENTAL) 

KA  =  THERMAL  CONDUCTIVITY  OF  THE  FLUID  (AIR) 

KP  =  THERMAL  CONDUCTIVITY  OF  THE  SPINE  MATERIAL 

M  =  A  HEAT  TRANSFER  PROPERTY 

NU  =  NUSSELT  NUMBER  (FROM  EXPERIMENTAL  DATA) 

NUC  =  NUSSELT  NUMBER  (FROM  CHURCHILL  AND  CHU  CORRELA¬ 

TION) 

NUAVG  =  AVERAGE  NUSSELT  NUMBER  (FROM  EXPERIMENTAL  DATA) 
NUT  =  TOTAL  OF  NUSSELT  NUMBERS  FROM  16  SPINES  ( 1  RUN ) 

PHI  =  ROLL  ANGLE  (DEGREES) 

PR  =  PRANDTL  NUMBER  (EXPERIMENTAL) 

R  =  RATIO  OF  TEMPERATURE  EXCESSES  (  =  THB/THO) 

RA  =  RAYLEIGH  NUMBER  (EXPERIMENTAL) 

RAT  =  TOTAL  RAYLEIGH  NUMBER  FOR  16  SPINES  (1  RUN) 

RAV,  =  AVERAGE  RAYLEIGH  NUMBER  (EXPERIMENTAL) 

TB  =  BASE  TEMPERATURE  (DEGREES  F)  UNCORRECTED 
TO  =  TIP  TEMPERATURE  (DEGREES  F)  UNCORRECTED 
TBF  =  CORRECTED  BASE  TEMP.  (DEGREES  F) 

TOF  =  CORRECTED  TIP  TEMP.  (DEGREES  F) 

TB ( I )  =  CORRECTED  BASE  TEMPERATURE  (KELVIN) 

TO ( I )  =  CORRECTED  TIP  TEMPERATURE  (KELVIN) 

THB  =  TEMPERATURE  EXCESS  AT  THE  BASE  OF  SPINE 

THO  =  TEMPERATURE  EXCESS  AT  THE  TIP  OF  SPINE 
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TF  =•  TEMPERATURE  OF  FILM  BETWEEN  SPINE  AND  FLUID 
TS  =  TEMPERATURE  OF  SURROUNDING  FLUID  (AIR) 

*###*###** ■*##-»#*DECLARATION  of  VARIABLES**************** 

REAL  D , KP , B , TS , R , M , H , PHI , NU , THB , THO , RAT , RAV , NUT , NUAVG , 
REAL  A , C , E , G , TF , VA , ALFA , KA , GR , PR , RA , NUC 
REAL  TB(25) ,T0(25),CB(25), C0(25) 

INTEGER  I 

****#*##*****#];NITIALIZATION  of  VARIABLES*************** 

D  =  0.0127 
B  =  0.0508 
KP  =  120 
G  =  9.807 
TS  =  87.53 
PHI  =  45 
NUT  =  0 
RAT  =  0 

CALCULATE  TS  IN  KELVIN 

TS  =  (TS  -  32 ) *5/9  +  273.15 

PRINT  HEADINGS 
WRITE (6,1 00 )  TS,PHI 
WRITE(6,1 50) 

READ  IN  AND  CALCULATE  CORRECTED  TEMPERATURES,  DEGREES  F 
DO  50  I  =  1,16 

READ (5i200)TB(l) ,T0(I) ,CB(l) ,C0(l) 

TBF  =  TB ( I )  +  CB ( I ) 

TOF  =  TO ( I )  +  CO (I) 

CONVERT  CORRECTED  BASE  &  TIP  TEMP.  TO  KELVIN 

TB ( I )  =  (TBF  -  32)*5/9  +273.15 
TO ( I )  =  (TOF  -  32 ) *5/9  +273.15 

CALCULATE  THETA  BASE  AND  THETA  ZERO  (TEMP.  EXCESSES) 

THB  =  TB(I )  -  TS 
THO  =  TO ( I )  -  TS 

CALCULATE  RATIO  OF  THETA  B  TO  THETA  0  (=cosh(MB)) 

R  =  THB/THO 
CALCULATE  M 

M  =  (ALOG (R  +  SQRT( (R**2)-1 ) ) )/B 
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C  CALCULATE  HEAT  TRANSFER  COEFFICIENT  H 
C 

H  =  (M**2)*KP*D/4 
C 

C  CALCULATE  FILM  TEMPERATURE 
C 

TF  =  (1 ./4. )*(TB(I)  +  TO ( I ) )  +  (1./2.)*TS 
C 

C  CALCULATE  PROPERTIES  OF  AIR  AT  TF 
C 

KA  =  0.004033  +  (7.30E-5)*TF 

VA  =  (-1.22567  +  (9.350E-3)*TF)#(1 .OE-5) 

ALFA  =  (-1.82333  +  ( 0. 01 350*TF )  )*(1 .OE-5) 

C 

C  CALCULATE  NUSSELT  NUMBER  BASED  ON  EXPERIMENTAL  H 
C 

NU  =  H*D/KA 
NUT  =  NUT  +  NU 
C 

C  CALCULATE  GRASHOF,  PRANDTL,  RAYLEIGH,  AND  NUSSELT  NUM- 
C  BERS  BASED  ON  CHURCHILL  &  CHU  CORRELATION 
C 

PR  =  VA/ALFA 

GR  =  (G*(1 .0/TF)*(TF-TS)*D**3)/VA**2 
RA  =  GR  *  PR 
A  =  0.387*RA**(1 .0/6.0) 

C  =  (1.0  +  (0.559/PR)**(9. 0/16. 0))**(8. 0/27.0) 

E  =  A/C 

NUC  =  (0.60  +  E ) **2 
RAT  =  RAT  +  RA 
C  PRINT  OUTPUT 
C 

WRITE(6,250)I,TB(I) ,TO(l) ,RA,NU,NUC,H 
50  CONTINUE 
C 

RAV  =  RAT/16.0 
NUAVG  =  NUT/16.0 
WRITE (6, 300 )RAV, NUAVG 
C 

STOP 

1  00  FORMAT ( 1 1 ' ,1 5X, » TS= ' , F6. 2, 3X , ' PHI= ' ,F4.1 ) 

1  50  FORMAT (17X, 'PIN' ,4X, ' TB ' , 7X, ' TO ' , 7X, ' RA ' , 8X, ' NU ' ,8X, 

§  *NUC» ,7X, *  H* ) 

200  FORMAT  (4  (F  7.3) ) 

250  FORMAT (1  8X, 12 , 2X, F6. 2, 3X, F6. 2 , 2X, F7 . 2, 2X, F8 . 5 , 2X, F8. 5 
#  3X,F6.3) 

300  FORMAT (1 6X, 'RAV  = ' , F7. 2, 2X, ' NUAVG  =',F6.3) 

C 
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